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Abstract 

A global statistical analysis of the pre-SNO solar neutrino data including the rates 
and recoil electron spectrum in SuperKamiokande is made assuming the solar neutrino 
deficit to be resolved by the interaction of the neutrino magnetic moment with the 
solar magnetic field. Given the general characteristics of the field profiles known to 
lead to the best event rate predictions, several specific choices of profiles are assumed 
and global fits investigated. Whereas previous studies revealed an excellent quality of 
the rate fits within the magnetic moment solution to the solar neutrino problem, the 
global fits are found to be only of reasonable quality and comparable to the global 
oscillation fits. This is related to the difficulty of any theoretical model to provide at 
present an adjustable prediction to the spectral data. 
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Whereas the apparent anticorrelation of the neutrino event rate with sunspot activity 
claimed long ago by the Homestake collaboration |I[ remained unconfirmed by other exper- 
iments p|, 0, 0] and theoretical analyses f^, the magnetic moment solution to the solar 
neutrino problem is at present an important possibility to be explored in the quest for an 
explanation of the solar neutrino deficit. This is the idea, originally proposed by Cisneros 
P], and later revived by Voloshin, Vysotsky and Okun 0, that a large magnetic moment 
of the neutrino may interact with the magnetic field of the sun, converting weakly active to 
sterile neutrinos. It now appears in fact that this deficit is energy dependent, in the sense 
that neutrinos of different energies are suppressed differently. In order to provide an energy 
dependent deficit, the conversion mechanism from active to nonactive neutrinos must be 
resonant, with the location of the critical density being determined by the neutrino energy. 
Thus was developed the idea of the resonant spin flavour precession (RSFP) proposed in 
1988 [Q. It involves the simultaneous flip of both chirality and flavour consisting basically 
in the assumption that the neutrino conversion due to magnetic moment and magnetic fleld 
takes place through a resonance inside matter in much the same way as matter oscillations 
. A sunspot activity related event rate in a particular experiment would hence imply that 
most of the neutrinos with energies relevant to that experiment have their resonances in the 
sunspot range. However, the depth of sunspots is unknown (they may not extend deeper 
than a few hundred kilometers) and the observed fleld intensity is too small in sunspots to 
allow for a signiflcant conversion. The anticorrelation argument has therefore lost its appeal 
for several years now. 

Despite the absence of the anticorrelation argument, several main reasons may be in- 
voked to motivate RSFP and investigating its consequences for solar neutrinos. In fact both 
RSFP and all oscillation scenarios indicate a drop in the survival probability from the low 
energy (pp) to the intermediate energy neutrino sector (''Be, CNO, pep) and a subsequent 
moderate rise as the energy increases further into the sector. The magnetic fleld proflles 
providing good flts to the event rates from solar neutrino experiments typically show the 
characteristic of a sharp rise in intensity at some point in the solar interior, followed by a 
progressive moderate decrease , |Tl| . This is in opposite correspondence with the energy 



dependence of the probability in the sense that the strongest fleld intensities correspond to 
the smallest survival probabilities. Hence RSFP offers a unique explanation for the general 
shape of the probability, which naturally appears as a consequence of the fleld proflle. On 
the other hand, from solar physics and helioseismology such a sharp rise and peak fleld 
intensity is expected to occur along the tachoclyne, the region extending from the upper 
radiative zone to the lower convective zone, where the gradient of the angular velocity of 
solar matter is different from zero [l^, Furthermore, it has become clear fl^, [|T0[| , 



TT| that RSFP provides event rate flts from the solar neutrino experiments that are far 



better than all oscillation ones JTSf, |]T6l. Finally, there are recent claims in the literature for 



evidence of a neutrino flux histogram containing two peaks, an indication of variability 
pointing towards a nonzero magnetic moment of the neutrino. 

The aim of this paper is to present a global statistical analysis of three of the fleld 
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profiles proposed in |jl0|], all obeying the general features described above, on the light of 



the updated standard solar model (BP 2000) and the most recent data [|I], ^, |^ . This 
global analysis includes the rates and SuperKamiokande spectrum only, as no day/night 
variations are expected for RSFP. Since the information on the SuperKamiokande total rate 
is already present in the flux of each spectral energy bin, this rate will not be included in 
the analysis, thus following the same attitude as in ref . [l^ . The SuperKamiokande energy 
spectrum used here corresponds to their new data with a range E^. = (5.0 — 20)MeV, 
19 data points and 1258 days exposure time [0. Moreover, all Gallium measurements 
(SAGE, Gallex and GNO) are combined in one single data point. For each profile the 
local minima of the global x^; defined as x^gi = X^rates + spectrum^ investigated in 
the ranges Am2i = (0 — 10^'^)eV'^ (mass squared difference between neutrino flavours) and 
Bq = {0 — 300)kG (value of the fleld at the peak). This Am^ range is related to neutrino 
resonances in the region from the upper radiation zone to the solar surface, that is, the whole 
region where a signiflcant magnetic fleld is expected. The upper bound of 300 kG for the 
peak fleld (Bq) at the base of the convection zone is suggested by most authors |jl2[, [l^. The 
most poorly known solar neutrino flux, namely hep neutrinos, is taken as a free parameter, 
whereas the ^B flux is flxed at its BP'2000 value The x^gi local minima occur in the 
range Amli = (0.7 — 2.0) x lO'^eV"^, with a peak fleld value from approximately 40kG up 

1.07-1.46 where the BP'98 hep flux ^ is taken 



The values of x^gi at these local minima are found to be close to each 



-11, 



and fhep = (j){hep) 

(analysis) /0(/iep)(Bp/98) 

as the reference one 

other and around 20 for 18 d.o.f.. The neutrino magnetic moment is taken at fii, = 10 ""fiB- 
This means that, since the order parameter is the product of /ii, by the magnetic fleld, for a 
value of the fleld at the peak near 300 kG, a suitable magnetic moment solution to the solar 
neutrino problem may exist at the remarkably low value /i^^ = 1.4 x 10~^^fiB- The results 
are presented in terms of 2a (95%CL) and 3cr (99.73%CL) contours on the Aml^, Bq plane. 



The proflles investigated are the following: 
Profile (1) 

B = , 



X < Xr 



R R X-Xr 
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where x is the fraction of the solar radius and xr = 0.70, xc = 0.85. 
Profile (2) 
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with xr = 0.65, xc = 0.75. 
Profile (3) 

B = Bi 
B = 



B = 2.16 X 10^ , X < 0.7105 

21 



^ /x-0.75\ 
~ V 0.04 ) 

Bo 

cosh30(a; - 0.7483) 



0.7105 <x< 0.7483 
, 0.7483 < a; < 1 



(7) 
(8) 

(9) 



with Bo = 0.9985i. 

The ratios of th 
recoil electron spectrum in SuperKamiokande is now 



The ratios of the RSFP to the SSM event rates Rq^ are defined as before |]T0| and the 



R 



th 



(10) 



for 19 energy bins (j=l,...19) 0. Here (piiE) is the SSM neutrino flux for component i 
{i = hep, ^B) and the factor fi is the ratio of the i-th flux relative to the SSM one. For 
i = hep the reference flux is taken to be the BP'98 one with /j as a free parameter, while 
for i =^B, for which (psB^^'^^ = 08^^^'^°°°, /i = 1 is assumed. The quantity /(E^, E^) is the 
energy resolution function ||21| of the detector in terms of the physical (E^) and the measured 
[Ee) electron energy [Eg = T + me)- The lower limit of E^ is the detector threshold energy 
= Eeth with Eeth= 5.0MeV) and the upper limit is evaluated from the maximum 
neutrino energy Em ||T0[ 



For the lower [10| and upper integration limits of the neutrino energy one has respectively 



En 



T' + \T'^ + 2meT' 



Em = IbMeV {i =^B) , Em = IS.SMeV {i = hep). (12) 



The weak differential cross sections appearing in equation (17) are given by 



daw Gp'^rrig 
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[{gv + gA)^ + {gv-gA)Hi 



T 
E 



gA 



TUeT 



(13) 



for z/eC scattering, with gy 



2sin^9w, gA = \- For z/^e and Ure scattering. 



^ = ^[(^- - 3a)' + + gA)' (l - I)' - i9v' 



gA 



rrieT 



(14) 



with gv = —k + 2sin^9w, gA 



The partial event rates for each neutrino component in each experiment were taken from 
and the solar neutrino spectra from Bahcall's homepage |^^. The contribution of the 
hep flux to both the Gallium and Homestake event rates was neglected. The analysis for 
the ratios of event rates and electron spectrum in SuperKamiokande was done following the 
standard procedure described in ||10[ f\. The validity of this procedure and alternative ones 
for solar neutrinos is discussed in refs. 12311 



The ratios of event rates to the SSM event rates and the recoil electron spectrum normal- 
ized to the SSM one, both denoted by -R*^ in the following, were calculated in the parameter 
ranges Am2i = (0 — 10~^)eV^^, Bq = {0 — 30) x lO'^G for all magnetic field profiles and 
inserted in the definitions for the rates and spectrum, 

xites= E (K - Rn^'n Kutot)Y^ljRf, ^ (15) 



31 

ii j2=i 



with (Ga=l, Cl=2) 



19 

j 1,32 = 1 



The quantities i?^^'P in eqs. (15), (16) are directly read from tables 1, 11 respectively and 
the total error matrices a'^{tot) are derived from the definitions given in [|1^], using and 
the error bars in tables I, II. In performing the fitting one has 21 experiments (2 rates and 
19 spectral data) and three free parameters, namely the mass squared difference between 
neutrino flavours Am^i, the peak field value Bq and the hep flux parametrized by fhep, hence 
18 degrees of freedom. With no correlations between the rates and spectral errors one has 
for the global x^, 

2 _ 2 I 2 ('\7\ 

X, gl X, rates ' X spectrum' \ ' ) 

For each of the three profiles, eqs.(l)-(9), the local minima of x% in the above parameter 
ranges were first investigated. Five such minima were found for the equilateral triangle, 
profile 1, all with x^^/lSd.o./. ~ 20 (see table III). For the second profile, seven local 
minima occur with a x^i/^^d.o.f. ranging from 20 to 24 (see table IV) and for the third one 
only two appear with xli/^8d.o.f. = 21 — 22 (see table V). Each of these minima is obtained 
for an hep flux from 1.07 to 1.46 of its BP'98 value as also shown. In view of the closeness 
of the values for all these fits with the possible exception of three of them, it is hard to 
tell which values one should expect the data to prefer. It is clear on the other hand that, 
since the value of the field at the peak may be as large as 300 kG and this analysis was done 
using fii, = 10^^^/iB, a value Bq = 0{4:OkG) (table III) means the existence of a possible 

^Here only the main definitions and differences are registered. For the calculational details we refer the 



reader to ref.| f 
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solution to the solar neutrino problem with /ij^ as low as 1.4 x 10~^^/iB in consistency with 
astrophysical bounds Moreover, the fits with the larger Bq values are disfavoured, as 
they would imply a value of the neutrino magnetic moment close to 10~^^ fiB, in conflict 
with most astrophysical bounds p5[] . 

The allowed physical regions and stability of these fits were also analysed. The results 
are presented in terms of 2a and 3a contours (95% and 99.73% CL respectively) in the 
Amli, Bq plane. These contours enclose the regions for which < Xmin + with, for 
18 d.o.f., Ax^ = 28.87 (95%CL) and Ax^ = 39.17 (99.73%CL) |2§. They are displayed in 
figs.l, 2, 3 for each of the profiles along with the local minima. The dashed lines denote 
2a and the full lines denote 3a contours. The local minima chosen for reference, that is, the 
ones for which the contours are defined, are numbers II (profile 1), VI (profile 2) and XIII 
(profile 3). In view of the closeness of the x^ values at the minima, their choice is irrelevant 
in terms of the difference among contours belonging to different minima: the discrepancies 
affecting the parameters Am\^, Bq are hardly perceptible on the scale of the plots. Finally, 
the contours were worked out for an hep flux {fhep) fixed to its value at the chosen minimum 
of x^- 

As an example, the predicted spectrum (eq.(lO)) for profile 3 is shown in fig. 4 at the 
global best fit (XIII) superimposed on the SuperKamiokande data. For this case x% = 21.11. 
The moderate rise occuring in the theoretical curve for > 12MeV is the effect of the hep 
neutrinos. It is clear from this figure the inherent difficulty in fitting any model prediction 
to the experimental data, a fact obviously also present in statistical analyses of oscillations. 

To conclude, the analysis of prospects for the magnetic moment solution to the solar 
neutrino problem reveals global fits of reasonable quality. Previous analyses |T^, ||11||,|]14 



unambiguously indicated from the part of the data on rates a preference for profiles with a 
steep rise across the bottom of the convection zone in which vicinity they reach a maximum, 
followed by a more moderate decrease up to the solar surface. This leads to an excellence of 
rate fits alone which is not shared by the oscillation rate fits. Interestingly enough, this class 
of profiles is the most consistent one with solar physics and helioseismology [jl2|, |13| . On the 



other hand the global analysis presented here for rates and spectral data shows a quality 
of fits comparable to the oscillation ones. This originates from the relative poor quality of 
the spectrum fits caused by the shape of a function for which no theoretical model is able 
to provide at present a properly adjustable prediction. 

In the present analysis only time averaged data were considered and a fitting was made to 
a time constant profile 'buried' in the solar interior. If, on the contrary, the active neutrino 
fiux turns out to be time dependent, a situation most likely to be interpreted through 
the magnetic moment solution with a time dependent interior field, the present approach 
is obviously inadequate. Averaging the event rates over time implies disregarding possible 
information in the data which otherwise is available if different periods of time are considered 



T7| . The robustness of such a procedure will greatly improve with the accumulation of more 
data. 
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Experiment 


Data 


Theory 


Data/Theory Reference 


Homestake 


2.56 ±0.16 ±0.15 


7.7±li 


0.332 ±0.05 


Ga 


74.7 ±5.13 


129±^ 


0.59 ±0.06 §,§ 


SuperKamiokande 


2.4 ±0.085 


5.15±J;° 


0.465 ±0.052 § 



Table I - Data from the solar neutrino experiments. Units are SNU for Homestake and 
Gallium and 10^cm~^s~^ for SuperKamiokande. The result for Gallium is the combined one 
from SAGE and Gallex±GNO. 



Energy bin (MeV) i?^ 



5 < Ee < 5.5 





.436 


± 


0, 


.046 


5.5 <Ee<6 





.438 


± 


0, 


.024 


6<Ee< 6.5 





.435 


± 


0, 


.019 


6.5 <E^<7 





.438 


± 


0, 


.015 


7 < < 7.5 





.463 


± 


0, 


.015 


7.5 < < 8 





.483 


± 


0, 


.016 


8<Ee< 8.5 





.465 


± 


0, 


.017 


8.5 <Ee<9 





.438 


± 


0, 


.017 


9<Ee< 9.5 





.450 


± 


0, 


.018 


9.5 < < 10 





.455 


± 


0, 


.019 


10<E^< 10.5 





.442 


± 


0, 


.021 


10.5 <Ee<ll 





.407 


± 





.022 


ll<Ee< 11.5 





.455 


± 


0, 


.026 


11.5 < Ee<12 





.423 


± 





.028 


12 < < 12.5 





.422 


± 


0, 


.033 


12.5 < < 13 





.481 


± 





.041 


13 < < 13.5 





.431 


± 


0, 


.047 


13.5 <Ee<U 





.603 


± 


0, 


.065 


U<Ee< 20 





.493 


± 


0, 


.049 



Table II - Spectral energy bins in SuperKamiokande (1258 days) and the corresponding 
ratio of the experimental to the SSM event rate [0]. 
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Fit 




Bo (G) 


fhep 


xyisd.o.f. 


I 


6.68 X 10-9 


4.1 X 10^ 


1.21 


21.32 


II 


6.75 X 10-^ 


6.73 X 10^ 


1.19 


21.49 


III 


8.26 X 10-9 


1.34 X 10^ 


1.36 


20.22 


IV 


7.5 X 10^9 


1.7 X 10^ 


1.19 


21.37 


V 


1.0 X 10-^ 


2.32 X 10^ 


1.46 


19.95 



Table III - Local minima of x^i ioi profile 1 (eqs.(l)-(3)). The reference for which the 
2(7 and 3a contours are determined and drawn in fig.l is fit II. 



Fit 




BoiG) 


fhep 


xyiSd.o.f. 


VI 


1.22 X 10-^ 


9.4 X 10^ 


1.24 


21.08 


VII 


1.21 X 10-^ 


1.22 X 10''' 


1.07 


24.29 


VIII 


1.28 X 10-^ 


1.7 X 10^ 


1.31 


20.33 


IX 


1.3 X 10-^ 


1.99 X 10^ 


1.09 


23.50 


X 


1.38 X 10-^ 


2.46 X 10^ 


1.43 


19.92 


XI 


1.41 X 10-8 


2.76 X 10^ 


1.11 


23.25 


XII 


1.56 X 10-8 


3.23 X 10^ 


1.45 


19.88 



Table IV - Local minima of x^i ioi profile 2 (eqs.(4)-(6)). The reference for which the 
2(7 and 3a contours are determined and drawn in fig. 2 is fit VI. 



Fit AmljeV^) Bp {G) hep xyiSd.o.f. 

XIII 1.36 X 10-8 1.04 X 10^ 1.26 21.24 

XIV 1.82 X 10-8 1.90 X 10^ 1.17 22.71 



Table V - Local minima of x^i for profile 3 (eqs.(7)-(9)). The reference for which the 2a 
and 3(7 contours are determined and drawn in fig. 3 is fit XIII. 
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Figure captions 

Fig. 1. Fits corresponding to the local minima of x^i for profile 1 (stars) and the 95% and 
99.73%CL contours (dashed and solid lines respectively) with respect to fit II. Since the 
values of x^i these minima lie close to each other (see table III), the contours may be 
regarded as actually enclosing the common 2a and 3a allowed areas on the scale of the plot 
for all fits. 

Fig. 2. Same as fig.l for profile 2. The reference fit is VI. 
Fig. 3. Same as fig.l for profile 3. The reference fit is XIII. 

Fig. 4. The theoretical prediction of the SuperKamiokande recoil electron spectrum (solid 
line) superimposed on the data set (1258 days) for profile 3 at the global fit XIII. Here 
the value of for the spectrum is xL = 21.11. 
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